Shot noise in magnetic tunnel junctions from first principles 
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We compute the shot noise in ballistic and disordered Fe|MgO|Fe tunnel junctions by a wave 
function-matching method. For tunnel barriers with <5 atomic layers we find a suppression of the 
Fano factor as a function of the magnetic configuration. For thicker MgO barriers the shot noise 
is suppressed up to a threshold bias indicating the onset of resonant tunneling. We find excellent 
agreement with recent experiments when interface disorder is taken into account 
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The statistics of electron transport in mesoscopic sys- 
tems has been subject to intensive research in the last 
decades leading to important and useful insights [H, Hj] ■ 
In a two-terminal conductor with a time-dependent cur- 
rent I(t), the simplest measure is the noise power P(w) = 
/^ o (AJ(0)A/(t))e^ t df, where AI(t) = I(t)-(I) denote 
the instantaneous fluctuation from the average current 
and (• • • ) a time and statistical average. The shot noise 
S is the zero frequency limit of the noise power when the 
applied voltage \eV\ is sufficiently larger than the thermal 
energy fcgT. The classical shot noise characterized by an 
uncorrelated Poissonian process is given by the Schottky 
formula S = 2e(I) Q. Shot noise contains information 
about the charge of the elementary excitations, entangle- 
ment, wave vs. particle nature of electron transport, and 
provides a diagnostic for open transport channels [j]. 

Magnetic tunnel junctions (MTJs) with MgO barriers 
0, IH have great potential for applications in magnetic 
random access memory elements and high-frequency gen- 
erators Band structure calculations of isomorphic 
Fc|MgO|Fe layered structures predicted a large drop in 
the electric resistance when the relative magnetization 
direction of the two ferromagnets switches from antipar- 
allel to parallel [H, 12|. The subsequently observed large 
tunnel magnetoresistance (TMR) [H, 0| can be explained 
in terms of the symmetry matching of only the majority- 
spin states in Fe with the Ai-band of MgO, which is the 
by far least evanescent in the gap. The tunneling ratio 
of the majority spin electrons is therefore relatively high 
while minority spin states arc efficiently filtered out by 
the MgO barrier. However, a quantitative first-principles 
description of transport in magnetic tunneling junctions 
is complicated by defects. The chemical composition of 
the interface strongly affects the TMR [13l4l5j and vari- 
ous interfacial defects have been identified to reduce the 
TMR The I-V curves alone cannot discrimi- 
nate between the possible different origins that reduce 
the TMR. 

According to conventional wisdom shot noise in tun- 



nel junctions is classical [l|, in agreement with earlier 
experiments ljl, 13] • Recent evidence that shot noise in 
MTJs is suppressed in the parallel configuration came 
as a big surprise [2lj . In order to resolve this issue we 
present parameter-free calculations of shot noise in mag- 
netic tunnel junctions. We compute sub-Poissonian shot 
noise for the parallel magnetic configuration and explain 
the results in terms of highly-transmitting resonant tun- 
neling between states localized at the interfaces between 
the ferromagnet and insulator. The agreement between 
first-principles theory and experiments plj is quantita- 
tive when disorder is taken into account. These results 
provide strong evidence of coherent transport and (addi- 
tional) proof for the very high quality of the MTJs used 
in that study. 

According to the scattering theory of transport a two- 
terminal conductor subjected to a sufficiently small bias 
voltage V leads to a time-averaged electric current 



and shot noise 



(I) = ^| e y|Tr(ttt) 



S= ^|eV|Tr(rtrt+t) 

h v ' 
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where t and r are the matrices of the transmission and re- 
flection coefficients in the space of the transport channels 
of the leads to the scattering region. These equations be- 
come more transparent by making use of the distribution 
function p (T) = £ n S (T - T n ) of the eigenvalues {T n } 
of the transmission matrix T = tH, where T n <G [0, 1] : 



(I) = j \eV\ J p(T)TdT 

s = ¥ |ey| / p(t)t(i - T)dT = 2eF{i) > 
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where F < 1 is the Fano factor. For a conventional tunnel 
junction transmissions are small and p (T) is substantial 
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only for T « 1. We then may disregard the ~T 2 term in 
the integrand of Eq. ((4]) and classical shot noise corre- 
sponding to F — > 1 is recovered. Clearly, a suppression 
of the shot noise that would correspond to a Fano fac- 
tor significantly smaller than unity requires that p (T) 
is significant at transmissions close to unity. Indeed, be- 
low we find such highly transmitting states in MT Js with 
sufficiently thin barriers. 

Previous theoretical treatments of the statistics of 
quantum transport have been limited to simple mod- 
els. While these can be sufficiently accurate for, e.g., 
structures defined on a two-dimensional electron gas, the 
details of the electronic structure are essential to under- 
stand (nearly) ballistic MTJs 



11, 12, £2|. This Letter 



reports the results of material-specific first principles cal 
culations of the statistics of transport in Fe|MgO|Fe mag 
netic tunnel junctions as a function of magnetic configu 
ration, voltage bias, and interface morphology and com 
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pare theory with experiments by Arakawa et al. 
in particular the suppression of the Fano factor for the 
parallel configuration. 
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FIG. 1: (Color online) Schematic diagram of a 
Fe|MgO|Fe(001) junction containing 5 MgO monolay- 
ers. The magnetization Mi of the left lead coincides with 
the z-axis; the magnetization M2 of the right lead lies in 
the zz-plane with angle 9. The red and blue grids in the 
scattering region denote the O and Mg atoms, respectively. 
We allow for interfacial disorder (oxygen vacancies, cyan grid) 
at the Fe|MgO interfaces. The applied bias eV — fiR — fiL, 
where p,R and fiL are the chemical potentials of the right and 
left leads, respectively. 

We consider an MTJ consisting of a MgO barrier and 
two semi-infinite iron leads as shown in Fig. [TJ The elec- 
tric current flows along the (001) crystal growth direc- 
tion. We incorporate the small lattice mismatch between 
the leads and the barrier by a 3% compression of the 
MgO lattice constant. The self-consistent calculations 
are carried out with the tight-binding linear muffin-tin 
orbital (TB-LMTO) [H surface Green function method 
24|. Disorder is treated using the layer CPA (coherent 
potential approximation) [25[. The atomic sphere (AS) 
potentials serve as input to the second step, in which 
the transmission matrix is calculated using a TB-LMTO 
implementation [26jj. Disorder is modeled by large lat- 
eral supercells, distributing the self-consistently calcu- 
lated CPA- AS potentials randomly laycr-for-layer in the 



appropriate concentrations in as many configurations as 
necessary. Depending on the defect concentration, most 
of our results are based on lateral supercells containing 
72 (two times 6x6) or 128 (two times 8 x 8) Fe atoms 
per monolayer . 

TMR = R (AP) /R (P) - 1, where R (AP) [R (P)} is 
the electric resistance for the antiparallel (AP) [paral- 
lel (P)] configuration of the lead magnetizations. For 5 
MgO monolayers (L) at low bias TMR = 3580% for spec- 
ular interfaces, which decreases drastically to 250% when 
5.56% oxygen vacancies (OV, the energetically most fa- 
vorable defect) are introduced at both interfaces. The 
TMR for the ideal junction is consistent with published 
calculations [ll[ , while that for disordered junctions is of 
the same order of magnitude as found in experiments [f| . 

Based on the calculated scattering matrix, we compute 
the Fano factor for various junction parameters. The re- 
sults are shown in Table [I] together with the TMR. For 
thick barriers, the Fano factors are very close to unity, 
implying full Poisson noise as expected. As the bar- 
rier gets thinner, the Fano factor of the parallel config- 
uration is increasingly suppressed. For a 5 MgO layers 
junction with 5.56% interfacial disorder, the Fano factor 
is F P = 0.87(4) and F A p = 0.98(1) for both configu- 
ration, close to the experimental values Fp = 0.91(2) 
and F AP = 0.98(1) for the same thickness [2l|. We 
can identify the majority (^) and minority (|.) spin con- 



tributions to be F- 



" 



0.96 



F 



HI) 



S 



too 



/ (2e(, 



and Fp 



0.72, where 
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In order to trace the origin of the shot noise suppres- 
sion, we plot the distribution functions of the transmis- 
sion matrix eigenvalues Pp/ap{T) in Fig. [5] involving 
7 x 10 6 eigenvalues over the whole Brillouin zone. For 
P, we identify a few high values of T n , which accord- 
ing to Eqs. (|3|4p affect shot noise S more strongly than 
conductance G: 0.3% of the eigenvalues are larger than 
0.05 but contribute about 39% to G but 89% to the in- 
tegrand proportional to T 2 in Eq. (J4|, which suppresses 
S. The integrands proportional to T and T 2 are shown 
for each eigenvalue interval in the histograms of Fig. [2] 
The dashed bars indicate a larger statistical error caused 
by the small number of eigenvalues at high T„. For AP, 
only very few T n fall into the region between 0.05 and 
0.1, the rest (99.95%) are all less than 0.05. 

In Table [J we can see that interfacial defects are nec- 
essary to explain the observed shot noise suppression 
211 ] . The OV concentration-dependent Fano factor for 



the 5MgO junction can be found in Table HU The statis- 
tics of our supercell calculation is found to be good by 
comparison with a Green function formalism in which 
impurity scattering is handled by the CPA 27J. Fur- 



thermore, even though an OV concentration around 5% 
suppresses the Fano factor, a further increase leads to a 
dramatical recovery of the full shot noise. The full shot 
noise observed in earlier experiments (l9L l20j is therefore 
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TABLE I: Barrier thickness dependence of the Fano factor in Fe|nMgO|Fe MTJs for the parallel (P) and antiparallel (AP) 
configuration. The results in square brackets are obtained for disordered junctions with 5.56% oxygen vacancies at the interfaces, 
where the error bar is given in round brackets. 



Fano factor 


3MgO 


4MgO 


5MgO 


7MgO 


P 


0.64[0.65(2)] 


0.91(0.69(4)] 


0.97(0.87(4)] 


1.00(0.99(1)] 


AP 


0.94(0.77(2)] 


1.00[0.94(l)j 


1.00(0.98(1)] 


1.00[0.99(1)] 


TMR 


1320% [165%] 


2400% [288%] 


3580% [250%] 


5600% [107%] 



TABLE II: OV concentration dependence of the Fano factor in 5MgO MTJs for P and AP configurations. The impurity 
concentrations are obtained by different numbers of OVs at the interfaces per lateral unit cell. 2 OVs in the 6x6 supercell 
correspond to 5.56%, 3 or 4 OVs in an 8 x 8 supercell correspond to 4.69% and 6.25% respectively. The conductances are given 
m units of 10 _5 e 2 //i per Fe atomic interfacial area, where the statistical error bar is given in round brackets. The conductances 
in square brackets are obtained by the Green function method [53] as a check of the statistics of the supercell method. 



Concentration 


P 


AP 




G(maj) 


G(min) 


Fano factor 


G 


Fano factor 





68.00[68.50] 


3.47(3.51] 


0.97 


1.95(1.95] 


1.00 


4.69% 


89(3) [94.5] 


41(7) [33.6] 


0.85(3) 


35(3) [32.8] 


0.98(1) 


5.56% 


80(4) [91.0] 


44(10) [29.2] 


0.87(4) 


36(4) [37.5] 


0.98(1) 


6.25% 


79(3) [85.5] 


25(4) [26.9] 


0.95(2) 


38(3) [41.0] 


0.98(1) 
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FIG. 2: Distribution of the transmission matrix eigenvalues 
p(T). The open circles are obtained for the AP configuration, 
for which 99.95% of eigenvalues are T n < 0.05. Solid trian- 
gles represent results for the P configuration. While most 
eigenvalues are still less than 0.05, the few high values prove 
the presence of resonant tunneling states. The red and blue 
histograms indicate the contribution of the integrands T and 
T 2 to each eigenvalue interval (for P). Dashed bars indicated 
increased statistical error due to the small number of T n 's 
approaching unity. 



consistent with higher disorder. 

In order to understand the sensitivity to the OVs, we 
plot the energy dependence of the conductance of 5MgO 
MTJs with different OVs concentrations in Fig. 03 In 
ballistic junctions the minority spin conductance for P 
and the AP conductance are strongly suppressed. Be- 
low the Fermi energy high transmissions are observed, 



FIG. 3: Energy-dependent conductance for P (minority spin) 
and AP configuration of a 5MgO MTJ for different impurity 
concentrations. Filled squares represent the ballistic junction, 
triangles 4.69% OV disorder, and circles 6.25% OV disorder. 
Red symbols stand for the minority spin conductance in the 
P and blue ones for the conductance in the AP configuration. 



however. For P these are caused by the overlap between 
interface states on both sides of the barrier. For AP, the 
interface states exist only on one side of the barrier, but 
since their symmetry is not orthogonal to the Ai states 
in the barrier and the majority spin states on the an- 
other side, the conductance is still high. Small amounts 
of OVs broaden and shift highly transmitting resonant 
channels toward the Fermi energy, thereby suppressing 
Fp (and the TMR). However, a further increase of the 
disorder destroys the resonant channels thereby recov- 
ering the full shot noise. The AP peak disappears and 
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FIG. 4: Angular-dependent Fano factor (triangles) and 
conductance (squares, per Fe atom at the interface) of 
5MgO MTJs with 5.56% OV disorder at the interfaces. 
When 8 changes from P to AP, the conductance decreases 
monotonously and Fano factor increases. 



FIG. 5: Bias voltage-dependent Fano factor of an antiparallel 
5MgO MTJ. I and S stand for current (per Fe atom at the 
interface) and shot noise integrated over the bias window, 
respectively. 



becomes a step structure near the Fermi energy. 

A 5% OV concentration appears to be close to the ex- 
periment [21| . since it explains both Fano factors and 



the TMR. Further information may be gained by the 
0-dependence of the Fano factor in Fig. @] for 5MgO 
junctions with 5.56% OVs. F (9) increases when moving 
from P to AP which can be understood by the arguments 
above. 

An important issue of the current- induced spin transfer 
torque (STT) in MTJs is its bias dependence. Recent ex- 
periments [28[ discovered a nonlinear increase of the STT 
and current at an applied bias of 0.2 V. This value if far 
below the MgO band gap as calculated in the local den- 
sity approximation (LDA) (22[. Since the LDA strongly 
underestimates band gaps, the observed threshold must 
have a different origin. A recent first principles analysis 
of the STT in Fe|MgO|Fe junctions explained the thresh- 
old in terms of resonant transmission channels in the AP 
configuration [29[. This hypotheses can be tested by the 
bias dependence of the shot noise. At finite bias V, the 
zero temperature current and shot noise read [l[ 



2 e eV 



I(V) 



p(T, E)TdT 



dE 



S(V) = j J p(T,E)T(l-T)dT 



(5) 



dE (6) 



Fig. [5] shows the integrated current and Fano factor as 
a function of applied bias for an Fe(f) |5MgO|Fe(4,) junc- 
tion. The Fano factor is unity for low bias but suddenly 
decreases with increasing bias at the threshold of the non- 
linear current characteristic, which for a clean junction 
is at about 0.8 V, consistent with the threshold bias in 
the STT. Small amounts of oxygen vacancies in MgO can 
lower this threshold bias to become closer to the experi- 



mental value [29[ . 

In conclusion, we compute sub-Poissonian shot noise 
for magnetic junctions with thin MgO barriers from first 
principles in good agreement with experiments. We in- 
terpret these results as strong evidence for resonant tun- 
neling states weakly broadened by disorder scattering. 
While it was known that MgO based tunneling junctions 
can be grown with high crystalline quality, we believe 
that the implied wave functions coherence over the tun- 
neling barrier is an important piece of information that 
lends credibility to the prediction of large thermal spin 
transfer torques j30[. This additional evidence of the su- 
perior electronic properties of MgO junctions should have 
ramifications for the application of this materials to other 
than ferromagnetic systems. 
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